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PREFACE 


The work described in this report was performed by the Systems 
Division of the Jet Propulsion Laboratory. 
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ABSTRACT 


The mathematical theory of the nutation of Mars is derived by 
classical rigid-body dynamics. The effect of nutation is to produce a 
26 -m maximum horizontal amplitude oscillation (at the surface of Mars) 
with a period of half a Martian year. This effect should be detectable 
in the Viking-Lander data. 


I. Introduction 

The nutation of the Earth results from the action of both the Moon 
and Sun and is a consequence of the small periodic parts of their gravi- 
tational couples (rhe nonperiodic parts produce the steady precession). 
Where the satellites of Mars are concerned, their motions are nearly 
in the equatorial plane, and thus their effect is negligible because of 
their small masses; consequently, it is sufficient to consider the action 
of the Sun alone. However, the eccentricity e of the Martian orbit is so 
much greater than that of the orbit of the Earth that the analytical 
expression for the solar couple on Mars mu, it be taken to higher order 

to achieve comparable accuracy. In what follows , all terms in the 

2 

nutation of Mars of order up to e will be found (e = 0. 0933 for Mars , 
as compared with 0.0167 for the Earth). 

The principal moments of inertia of Mars will be denoted by A, 

B, C, and chosen such that C > A, and that A = B. Numerically for 
Mars (C - A)/C = 1/192 (Ref. 1); the obliquity of the orbital plane in 
relation to the equator of the planet is 25° 12'; the sidereal rotation 
period 2tt/w is 24^^ 27 *^”^ 23®; and the orbital period 2Tr/n is 687 mean 
solar days. 
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II. Theory of Nutation 


This theory is given in adequate detail for present purposes in 
Chapter XX of Ref. 2, and provides the mathematical basis for this 
report. 

In Fig. 1, OX Y Z are fixed rectangular axes at the center of 
o o o 

Mars, with OX Y the plat.e of the orbit of the Sun about Mars. The 
o o 

principal moments of inertia (A, A, C) have axes OX, OY, OZ, 
respectively, as shown, and their positions relative to OX^Y^Z^ are 
defined by Eulerian angles 0, 4^, 4> as shown (Ref. 2 notation adopted). 
The direction of X will be define later. 


Z 

u 



Fig. 1. Mars Geometry 

It may readily be shown that the slightly elliptical charact3r of 
the equator of Mars, recently determined, would give rise to a fourth- 
order term in the potential U only af order 10”^ of that cn which the 
current dynamical calculation rests, and leads to differences in the 
values of surface measurements of lengths of a few tenths of a milli- 
meter at most. The fifth-order terms are less than 10~^^ of U 
itself. 
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Neglecting higher order effects, the Lagrange potential for 
determining the couples is shown (Ref. 2, p. 353) to be 

(I) 

n = mean motion of Mars about the Sun 
a = semimajor axis of orbit about the Sun 
p = Mars — Sun distance 

and the polar equation of the Martian orbit is 

a 1 + e cos V 

p “ , 2 

1 - e 


with V the true anomaly, while z in (1) is the perpendicular distance of the 
Sun above the Mars equatorial plane OXY. 

Under the assumption that 0, 0 and ij;, ijj , 0i|> are small com- 
pared withwijj and w0 the equations of motion leading to the nutation- 
terms (Ref. 2, p. 352) are, to sufficient accuracy. 


^ - 

-(C 0 ) sin 0)-^ 1^ 

(3) 

0 = 

(C 0 ) sin 0)”^ 

(4) 


In the steady precession of the planet, the angle increases uniformly 
with time, so that ij; = + pt (where 27r/p = period of precession), 

while 0 remains constant at a value 0^, say. It is the small differences 
from these 0^ in the motion of the point Z that constitute the 

nutation. 
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The function U needs to be expressed in terms of the time i and 
0, ifj^ and other constant angles associated with the orbital motion. 
But it is of interest that if U is developed in terms of the true 
anomaly v, then finite expressions emerge for 3TJ/3Q, dU/d^, and 
hence for 0, ij;, that would hold for all values of the eccentricity 
e < 1. To show this, it is readily found (Ref. 2, p. 358) that 

— = sin 0 sin (v + w + tjj) (5) 

P 

where w the angle subtended (at the center of Mars) from its vernal 

equinox^ to the Sun at perihelion. This is the same angle measured 

2 

between Mars* autumnal equinox and Mars' perihelion — centered 
at the Sun. 


Using the geometry and Ref. 3 for expressions for the orbital 
angles, we find approximately 

w = 250. 6804° + 0. 6343°T (6) 

where T = Julian centuries since 0^^ January 1, 1950. Here w includes 
a precession value (i. e. , takes into account some changes due to iji) 
and, further*, with i{J being so far indefinite in its zero point, it can 
be considered small. 

We then have 

3 

U = “!■ (C - A) n^ sin^ 0 sin^ (v H- w + v|j) (7) 

But in the elliptic motion 


2 dv , 2 2 

P = VI - e 


Direction from Mars to Sun at beginning of spring. 
Direction from Mars to Sun at beginning of Mars winter. 
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- z 

and hence d/dt = h p“ d/dv , so that Eqs. (3) and (4) can be transformed 
to involve differentiation with respect to v instead of t, thus 


dv 


-*.C to sin 0) 


-1 au 

no h 


(9) 


~= (C u sin 0)"^ 


8U 

94^ h 


( 10 ) 


2 

Since p is independent of 0 and ijj, the f*;:nction CJ' = (p /h)U can be used 
now instead of U, and we have 


U< =-|(C - A) n (1 


e ) 


-3/2 


sin G (1 + 3 cos v ) 




cos (2v 



Thus, the equations of motion take the form 


( 11 ) 


dv “ 


3 (C - A) 


_n 

to 


(1 - e 


2, -3/2 


cos 0 < 1 “ cos (Ev + 2w) T e cos v 


e cos (2\v + 3v) 


(2w + v)j- 


(12) 


and 


dO _ 3 (C - A) n ,, 2. -3/2 

_ . -j— ^ w - e ) 


sin 0 


sin (2w + Zv) +-^ e sin {2w -t- 3v ) 



e sin (2w + v 


>) 


( 13 ) 
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Since w increases with time (through ijj) only extremely slowly, it 
follows that insofar as (1) is adopted as adequate representation of 
the potential, the values of *P(v) and Q(u) obtained on integration of 
(12) and (13) are closed expressions valid for all values of e < 1. 

Thus it readily follows, omitting the very slow prccessional change 
of ijj with time, that 

^ cos 0 |v sin (2w + 2v) 

+ e sin V - sin (2w+ 3v) -■|■e sin (2w + v)| (14) 

and 

0 - 0^ = j ^ -~(1 - e^)’^/^ sin 0 cos (2w + 2v ) 

1 

+ •“6 COS (2w + 3v) +'2'® cos (2w+ t<)| (^5) 

In order to express 4i and 0 in terms of the time explicitly, it 
would be necessary to expand the trigonometrical terms in Eqs. (14) 
and (15) in terms of the time, which can readily be done by means of 
standard formulas. But for the purpose, series-expansions in powers 
of e are necessary and would hold only for sufficiently small values of 
the eccentricity. In fact, for Mar*; e = 0. 0933, and numerical 
results later obtained herein show that expansion as far as terms in 
e will give quite sufficient accuracy for any likely application, such 
as locations of positions on the Martian surface. Alternatively, the 
function U as given by Eq. (7), can be developed at once in terms of 
the time as a power series in e, and then use made of Eqs. (3) and (4). 
Both methods lead, of course, to precisely the same final results. 
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When the second procedure is adopted, since using M for the 
mean anomaly nt + c - S, 


5 2 

V = M t 2e sin M 4 ~ e sin EM 


it is found that 


U =-|-(C - A) sin^O 1 +|-e^ - cos (2M t 2w) 


9 2 1 7 

+ 3e cos M cos 2M + cos (M f 2\v) - cos (3M +2w) 


17 2 o ,v1 

— e cos (4M f Ev/Jj 


(16) 


Then forming Eqs, (3j and (4) from this expression for U, and 
integrating with resp jct to t, it is found that 


_ 3 (C - A)n ^ r/, . 3 2\ , 

= 2 — 0 ^1 + 2 ® 

/ 1 5 2 \ 1 

jsin(2M + 2w) + 3e sin M + -g- e sin (M + 2w) 

- e sin (3M + 2w)+ ~ sin 2M - sin (4M + 2\v) 


and 


(17) 


^ _ 3 (C - A) n . „ r/l 5 2\ 

G - 0^ ='T ' — — ' — sin 0 hr - ) 

o2 C to [_(2 4/ 


cos (2M + 2w) 


1 7 

- 2 ^ cos (M 4 2w) + ^ e cos (3M 4 2w) 


,11 2 

4 e cos (4M 


4 2w)j 


(18) 
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We now choose X^, so far an arbitrary direction, such that 
= 0 at the chosen epoch. Thus ip is always small and can be 
ignorec in the trigonometric arguments. 

The term proportional to nt in Eq. (17) represents the steadv 
3 

precession, and if we write i{j = + pt + nutation (in if^), the numerical 

values for the coefficients of the nutational terms in ijj and 0 are fou-d 
to be (in seconds of arc) as follows 

ijj = 4 ,^ + pt - 1'.’0715 sin (2M + 2w) + 0V6133 sin M 

+ 0'.'1022 sin (M -t 2w) - 0V2385 sin (3M + 2w) 

+ 0V0215 sin 2M - 0'.'0405 sin (4M + 2w) (19) 

0 = 0^+ 0'.'5042 cos (2M + 2w) - 0'.'0481 cos (M t 2w) 

+ 0'.'1089 cos (3M 4 2w)4 0'.'0185 cos (4M I- 2w) (20) 

It is seen that the largest coefficients are those of the terms in 

2M and thus have a period of half the Martian year. The greatest semi- 
amplitudes of sin 6 and 0 are seen to amount in sum to about 2 " in 
ijj at the equator, and 0'.'6 in 0. These would correspond, respectively, 
to displacements of about 26 m and 9 m of the surface relative to a 
steadily processing Mars. 

[Note added in proof, August 1979: Analysis of the Viking-Lander 
data has confirmed oscillations of the foregoing order of magnitude. 

The relevant data and analysis will be published in due course.] 


^The value of precession for this development^ 7'.’415/yr, follows 
from (C - A)/C = 0. 00524 (Ref, 1), 
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